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Abstract 

This work addresses the question of the relation between strike-point splitting and magnetic stochasticity 
at the edge of a poloidally diverted tokamak in the presence of externally imposed magnetic perturbations. 
More specifically, ad-hoc helical current sheets are introduced in order to mimic a hypothetical screening of 
the external resonant magnetic perturbations by the plasma. These current sheets, which suppress magnetic 
islands, are found to reduce the amount of splitting expected at the target, which suggests that screening 
effects should be observable experimentally. Multiple screening current sheets reinforce each other, i.e. less 
current relative to the case of only one current sheet is required to screen the perturbation. 
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1. Introduction 



An axisymmetric, single null, poloidally diverted tokamak has two strike-lines where the plasma hits 
the divertor targets: one on the high field side (HFS) and one on the low field side (LFS). In the presence 
of a non-axisymmetric magnetic perturbation, these strike- lines are replaced by spiralling patterns [![. If 
an experimental profile (e.g. D Q or infrared [IR]) is taken along the radial direction at a given toroidal 
location, the strike-points are then observed to split. Such a splitting is for instance commonly observed 
during locked modes Q. In the presence of external resonant magnetic perturbations (RMPs) a splitting 
may also be expected. On DIII-D, the splitting is observed during Edge Localized Mode (ELM)-suppressed 
discharges using n = 3 perturbations from the I-coils [3, 4]. It is important to notice that the splitting is seen 
much more clearly on particle flux (D Q ) profiles than on heat flux (IR) profiles, at least in low collisionality 
experiments Q (at high collisionality, the splitting is however clearly observed on the heat flux Q). Recently, 
DIII-D also reported on splitting observations (both on heat and particle fluxes) in L-mode plasmas [S] . JET 
(using the Error Field Correction Coils (EFCCs)) and MAST (using n = 3 perturbations from the ELM 
control coils (|) find consistent effects on the heat flux profiles: the splitting is observed in L-mode but not 
in H-mode 0]7 

In the field of ELM control by RMPs from perturbation coils, one major question is to know whether 
the RMPs stochastize the magnetic field at the edge of the plasma, as assumed by the vacuum modelling. 
Studies based on the vacuum field assumption [§] have led to a design criterion for the considered ITER 
ELM control coils However, two important elements cast doubt on the stochastization of the magnetic 
field. The first one is the absence of a degradation of the electron temperature gradient in the edge transport 
barrier, which would be expected in the presence of a stochastic field [10]. The second one is the strong 
rotational screening effect [THEl found in simulations of the DIII-D ELM suppression experimen ts I 13L Il4| . 
On the other hand the RMPs can also become amplified by the interaction with MHD modes [15| — an 
effect which works against rotational screening [lq ]. 

In this paper, we analyse the possible consequences of the rotational screening (without taking into 
account the possibility of amplification by MHD modes) on the strike-point splitting in order to assess 
whether screening effects may explain the absence of a clear splitting of the heat flux profiles in some 
experiments, in particular in the DIII-D, JET and MAST H-modc discharges referred to above. 



2. Modelling and theoretical understanding of the strike-point splitting 

Under a non-axisymmetric perturbation the magnetic separatrix splits into two surfaces: the stable 
and unstable manifolds of the X-point. The stable (resp. unstable) manifold is the set of field lines that 
asymptotically approach the X-point when followed in the direction of (resp. opposite to) the magnetic field. 
The manifolds are of interest to experiments because they delimit the first passage through the wall of field 
lines arriving from the plasma core. Their intersections with the divertor plates thus define areas (divertor 
footprints) where high heat and particle fluxes are carried from the plasma core along the field lines (l7j . 
Those areas take typically the form of spirals of high temperature and particle recycling around the original 
(unperturbed) divertor strike point. 

The divertor footprints can be visualized by plotting a map of the connection length on the divertor 
plates (a laminar plot) [18j . The connection length is the distance (measured as the number of toroidal 
turns) needed to reach the wall again by following a field line starting at a given position. Field lines with 
large connection lengths remain in the plasma for many turns and carry high fluxes from the hot plasma 
core. The extent of the footprint can be approximated analytically using the Melnikov function II QB whose 
maximum is the difference of ip between the unperturbed strike point and the tip of the footprint [20| . When 
the perturbation has one dominant toroidal mode further simplification is possible and the difference of ip 
can be expressed using a single number — the one- mode Melnikov integral M n [2fj| . 

An example of the laminar plot and the stable manifold is shown on Figure[T](left plot) for an equilibrium 
predicted for the COMPASS tokamak [2l[ in the case of a magnetic field of 1.2 T, low triangularity, single- 
null (SND) geometry and heating by one co-injected neutral beam [22]. The n = 2 perturbation is imposed 
by the existing perturbation coils whose description can be found in [23| . 
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3. Physics-motivated method for taking screening currents into account 



3.1. Coordinate system and resonant field components 

We use an (s, ip, 6*) system of equilibrium coordinates, where s = — V'axis) / (V'sep — V'axis)! 1 ^ 2 (with 
ip the poloidal magnetic flux), tp the geometric toroidal angle and 9* the corresponding straight field line 
poloidal angle (i.e. such that \dip/dd*\ = const. = q along a field line, with q the safety factor). 

Magnetic islands are known to arise from the component of the magnetic perturbation which is perpen- 
dicular to the equilibrium flux surfaces. We characterize the latter by the quantity b 1 = B ■ Vs/B ■ Vip. It 
can be shown that its Fourier components b mn are directly related to the half-width of the magnetic islands 
and moreover they are related to the Melnikov integral |24| , more precisely proportional to a function which 
is a generalization of the one-mode Melnikov integral M n (the Melnikov-like function) [2(| or equivalently 
the Poincare integral (25|. At the same time the radial extent of the invariant manifolds (estimated by 
the Melnikov integral) gives the lower bound of the stochastic layer width because the intersections of the 
invariant manifolds (the homoclinic tangle) create themselves a thin stochastic layer. The actual stochastic 
layer can be much wider than this lower bound because it is formed also by the overlap of the magnetic 
islands [Hj]. 

3.2. Model of the screening currents 

Without loss of generality, we consider one toroidal mode n of the screening currents in the plasma. A 
generic current can be represented a Fourier sum of these modes which are independent due to the toroidal 
symmetry. 

Screening currents are modelled under the following assumptions: 

1. They are radially localized on infinitesimally thin layers around the resonant surfaces: j = X) 9 eS — 
Sq,n)jq,n, where S is the Dirac delta function and S is the set of rational values of the safety factor 
which define the screening surfaces: q — m/n for integer m and the given toroidal mode number n. 
The corresponding values of the radial coordinate s are noted as s 9)n . 

2. They are parallel to the equilibrium field lines: j q:n = j q , n / B oq ■ B eq . 

3. They are divergence-free: V • f — 0, which implies that a q ^ n = j q , n / B eq is constant on a field line. 
The first assumption corresponds to the fact that the current density is generally localized in a thin layer 



around the resonant surface [ll|, |12|, |26| . The second assumption follows from the fact that the screening 
currents are induced to oppose the radial perturbation and to create a radial screening field perpendicular 
to the field lines a parallel current is needed. The third assumption expresses quasi-neutrality. 
The angular dependence of a q ^ n has the form of one Fourier mode: 

a q , n (0*,<p) = K (/3 m „ exp [i (m0* + nip)}) (I) 

with m = qn, where the toroidal dependence is the consequence of working with one toroidal mode, while 
the poloidal dependence follows from the third requirement: a q ^ n constant on field lines. /3 mn is a complex 
quantity containing both the amplitude and phase of a q ^ n . Thus, the screening current density can be 
expressed as a linear combination of basis currents J mn ,o 

J = ^ ImnJmnfi , (2) 

\m;m/n£S J 

with coefficients I mn = (3 m nB TC f/j TC f and the basis currents 

J mn ,o (s, 0*,<p) = -^r~S (s - s qtTl ) exp [i (m9* + nip)} ■ B cq . (3) 

£>ref 

-Bref/jref is an arbitrary value expressing the choice of basis current amplitudes (and thus the normalization 
of I m n) relative to the magnetic field strength. 
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3.3. Coupling matrix 

We calculate the field B mn fl created by the base resonant current j m n,o- The numerical method ap- 
proximates jmn.o by discrete helical current filaments on the screening surface and calculates their vector 
potential on a flux surface aligned mesh. It avoids the mesh points which are too close to the screening 
surface and uses instead a cubic spline interpolation in order that the discrete approximation of currents 
not cause an error. B rnn _o is then obtained as the curl of the interpolated vector potential, which gives all 
its components and automatically satisfies the condition of zero divergence. The corresponding b 1 compo- 
nent, denoted &J' mra , is then Fourier transformed at each resonant surface q — m' /n' in order to obtain the 
resonant components - Thus, fr^™? designates the resonant part, on the q = m'/n' surface, of the 
b 1 created by a resonant current j m n,o located at the q — m/n surface. Due to the toroidal symmetry of the 
field equations in 6^™? we nave n = n', otherwise ^^™"o = 0- 

The plasma response field corresponding to the total current j as given by Eq. ([2]) is -Bpiasma = 
^ (E m ;m/nes ImnB m „,o) whose resonant b 1 components are 6^/ n , plasma = E m ^""o' 7 ™"/ 2 and ^-m'-n.plas! 

Tlmi (^m™"o ' Ann J /2- The RHS is the product of the matrix fo^™" with subscripts m and m', which we 
call the coupling matrix, by the current vector I mn . 



3.4- Calculation of the screened field by inversion of the coupling matrix 

To determine the coefficients I mn one needs an assumption about the character of the plasma response, 
i.e. if it amplifies or screens the perturbation and by what amount. In the following we assume an efficient 
screening which completely eliminates magnetic islands at the rational surfaces in the pedestal region, i.e. 
the resonant Fourier components b\ m of the total magnetic field are zero. 

The procedure to obtain the screened field begins with the calculation of the coupling matrix, for a given 
choice of the set of screening surfaces S. Independently, the vacuum RMP spectrum 6j n ,„, vac is calculated 
from the coil geometry. The screening current distribution Zj^ ccn is obtained by solving J2 m ^mW^o ' 7 mn C ° n = 
~2b] n , n vac , i.e. by inverting the coupling matrix. The full, screened field -Bfuii is then obtained as -Bfuii = 
-Svac + K (j2 m Imn^Bmn.jjj ■ It is easy to verify that its resonant Fourier components fr^^full on rational 

surfaces with q = m/n £ S satisfy the property b\ nn full = up to the error introduced by the numerical 
method. 

If the set S consists of only one screening surface, the coupling matrix is trivial, with one element. The 
screening currents for each surface alone are given by using only the diagonal terms of the coupling matrix: 
Imndiag = _ ^Jnn,vac/^mnO' If more screening surfaces are considered, the currents I™ on each one 
required to cancel the resonant components may be different than for one surface alone due to off-diagonal 
components of the coupling matrix. This effect can be quantified by the ratio between i^ con calculated 
using the full matrix and I™" . If for example the ratio is lower than 1, it means that the perturbation 
field from different screening surfaces reinforce each other and smaller currents are required for screening 
than if screening were due only to one surface. This is the case for all the examples described in the next 
subsection, where we present actual values of this ratio. This is one of the geometry effects neglected in 



cylindrical models such as [13, [1J, [27[ where the coupling matrix is always diagonal. 



3.5. Numerical examples 

We chose two cases to illustrate the method and to show the effect of screening on footprints: the 
COMPASS case described above and an equilibrium from an ELM control experiment on JET with n = 2 
perturbation of the EFCCs (shot #79729 at 19.38s). For both cases we first calculate the screening field 
needed to cancel the perturbation on a single surface (q — 4/2 for COMPASS, q = 5/2 for JET). In those 
cases the coupling matrix is trivial, with one element. Then we calculate the screening field choosing four 
screening surfaces with q = 4/2, . . . , 7/2 for COMPASS, q = 5/2, . . . , 8/2 for JET. The n = 2 mode of the 
screening field 6g Croon as a function of 9* on the innermost screening surface is shown on Figure for the 
four cases (each equilibrium with one and four screening surfaces). The field is given on innermost surface 
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Figure 1: Laminar plot of the connection length (as a number of toroidal turns) on the divertor of COMPASS near the HFS 
strike point with a vacuum (left) and screened (right) perturbation field. White line: The stable manifold of the vacuum 
perturbation field. Black line: the stable manifold of the screened perturbation field. 

in order to show how the field changes depending on the choice of screening currents while maintaining the 
same resonant component. The screening field of one screening surface is distributed all over the resonant 
surface, due to the helical structure of the screening current, while for four screening surfaces the screening 
field is mostly localized on the LFS. Both COMPASS and JET show this effect. It shall be noted that 
the vacuum field is also localized at the LFS because of the position of the coils in both tokamaks, so the 
screening field of four currents is more similar to the vacuum field than the field of a single current. The 

ratio /S =2 A™ =2 ,dia g fOT JET is °- 63 ' fOT COMPASS 1^=2 A™ =2 ,dia g is 0.72, with similar and 
generally decreasing values for other surfaces. 
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Figure 2: Poloidal dependence of the real part of the n = 2 component of the relative perturbation field ft 1 at the q = 5/2 
resonant surface for JET, shot #79729 and q = 4/2 for COMPASS. Fields of one screening current (full lines) and four screening 
currents (dashed lines) are shown. 

The numerical error in determining the screening current amplitudes was estimated by increasing four 
times the resolution of the discretization of the screening currents (from 4096 to 16384). The largest change 
was observed for the outermost surface (6.3% for COMPASS, 8.3% for JET). The error decreases for the 
inner surfaces down to 0.1% for COMPASS, 0.2% for JET at the innermost surface. 
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4. Impact of screening currents on the splitting 



For the &J' mn field created by the resonant current j* m n,o we may compute the one-mode Melnikov integral 
at the separatrix which we will note M™ n . The total Melnikov integral which estimates the splitting from 
all the screening currents and the vacuum field is 



M n = M„, vac + ^ 



rscrccn ]\ / rm 
i mn 1V1 



(4) 



The ratio of the footprint extent measured in terms of tp for the screened vs. the vacuum field is given by 



Mr, 



M n 



Table 0] lists this value for the COMPASS and JET cases as a function of the choice of 

screening currents, and Figure [3] shows the same data. 

In both cases a significant reduction of the footprints is predicted by the Melnikov integral when four 
screening currents are considered. To confirm this we plotted the stable manifold and a laminar plot around 
the inner strike point for COMPASS for the screened field and the stable manifold for the vacuum field for 
comparison (Figure [TJ right plot). The stable manifold forms the boundary of the footprints as expected and 
indeed shows a clear reduction in comparison to the vacuum stable manifold. The difference of tf> (normalized 
to the poloidal flux at the separatrix) between the footprint tip and the base is 6.8 x 10~ 3 , while the Melnikov 
integral method predicts 5.3 x 10 -3 . The inaccuracy of the Melnikov method is significantly lower (relative 
error below 1/10) for the other cases in the table which do not include the outermost (m = 7) screening 
surface. Similar result was found for the JET case (laminar plots are shown in |7J): the actual difference is 
3.4 x 1CP 3 , while the Melnikov integral method predicts 3.6 x 10~ 3 . The outermost (m = 8) surface can also 
lead to a large error in the Melnikov integral estimation, especially with the m = 8 surface alone, where the 
Melnikov integral predicts reduction of footprints by a factor of 0.89, while the actual reduction is 0.56. 



COMPASS 


m — A 


m = 4, 5 


m = 4, 5, 6 


m — 4, 5, 6, 7 


|M 2 |/|M 2)Vac | 


0.73 


0.54 


0.45 


0.27 


JET 


m = 5 


m = 5, 6 


m = 5,6,7 


m = 5, 6, 7, 8 


|M 2 |/|M 2)Vac | 


0.69 


0.50 


0.30 


0.26 



Table 1: One mode Melnikov integral of the screened perturbation normalized to the Melnikov integral of the vacuum pertur- 
bation, for different choices of the screening currents, given by the resonant poloidal mode numbers m. |-M2,uacj is 2.00 X 10~ 2 
for COMPASS, 1.38 X 10~ 2 for JET, relative to the poloidal flux at the separatrix. 



5. Discussion and conclusions 

We developed a model of the plasma response currents on resonant surfaces and the resulting field, based 
on the realistic geometry of poloidally diverted tokamak plasmas and thus appropriate for the region near 
the separatrix, which is crucial for the ELM mitigation by external perturbations and also for the impact 
of perturbations on the divertor strike points (strike point splitting). To compute the screening currents we 
used the assumption of complete screening of resonant modes in the edge region, because we do not simulate 
the plasma response self-consistently. This is justified by earlier results indicating that the resonant modes 
of the perturbation will be suppressed by strong pressure gradient in the pedestal region [13l . [14| . As we do 
not directly couple our model with these results yet, we performed a scan of different possible combinations 
of screening currents. In future work we plan to use the results of self-consistent MHD models to determine 
the set of screening surfaces and the screening factors, which will also allow for incomplete screening. The 
reason for this approach is that those MHD models use a cylindrical geometry in which the effect on the 
divertor strike points can't be represented, and our model fills this gap. While there are MHD codes (e.g. 
NIMROD [lj|, M3D (2|| and JOREK [29[) using a realistic geometry which thus can model the strike points 
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Figure 3: The relative Melnikov integral from Table 1 as a function of the number of screening surfaces: from one at the q = 5/2 
resonant surface for JET, shot #79729 and q = 4/2 for COMPASS, up to four of them. 



themselves (as shown for JOREK in [29j) without needing our model, they have limitations in the physics 
included (lack of realistic resistivity) which justifies the interest of our approach. 

The resulting screened field was used to model magnetic footprints on the divertor by tracing the field 
lines. For two example cases (single-null equilibria of COMPASS and JET with n = 2 perturbations) we 
have shown notable differences in comparison with the vacuum field. The screening significantly reduces 
the spiralling patterns of field lines coming from inside the plasma. The spirals are shortened along their 
axis, the position of the axis is not affected. Reducing the coil current in a vacuum model has a similar 
impact. Comparison with experimental observation of strike point splitting could be thus used for validating 
the starting assumption about the screening. Coupling with MHD models will also enable us to do scans 
of the dependence of strike point splitting on the collisionality (and thus resistivity) and rotation, which 
are both important parameters for the screening effect, and compare the results with experiments. It 
can be expected that higher collisionality and slower rotation will reduce the screening and enhance the 
footprints. The discussion of experimental results is however more complicated than a simple comparison 
with the predicted magnetic footprints because of the deformation of flux surfaces beyond the separatrix. 
The particles in the scrape-off layer would follow those distorted surfaces, and the screening may change 
their distortion differently than in the case of the separatrix. This shows that it is necessary to distinguish 
between particle and heat flux and to evaluate the impact of screening on flux surface distortion in the SOL 
in addition to the separatrix. Indeed, as mentioned in the Introduction, in the DIII-D low collisionality 
experiments the splitting of particle flux and heat flux are different. Transport modelling using the vacuum 
perturbation failed to explain this observation [30j | and it will be thus interesting to repeat it for the screened 
field as calculated by our model to see if the screening can provide an explanation. 

The reduction of footprints can be efficiently estimated by the Melnikov integral method, but it is 
sometimes inaccurate when a screening current very close to the separatrix is included, probably because of 
the strong variation of the screening field while the Melnikov method (a first-order perturbation method) 
needs the perturbation field to be slowly varying in the vicinity of the separatrix. 

The reduction of footprints can be qualitatively understood from the fact that the screening field is mostly 
localized at the LFS, just as the field of the coils. It can be shown that for a LFS-localized perturbation field 
the Melnikov integral which estimates the splitting is linked to the values of the Melnikov-like function on 
the resonant surfaces which is proportional to the resonant modes of the perturbation 20] . Thus eliminating 
the resonant modes shall also mostly eliminate the Melnikov integral and splitting. It shall be noted that for 
a single screening surface the field is not at all localized so this reasoning does not apply. Indeed, we have 
seen that a single screening surface reduces the Melnikov integral only weakly. Multiple screening surfaces 
also require less current to screen the resonant component of the field than a single screening surface — we 
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may say that the currents reinforce each other. This effect is due to the realistic geometry used here and 
can't be represented in the cylindrical approximation used in some models for numerical simplicity. 

Our model is designed to represent one significant feature of the plasma response to the perturbations 
— the surface screening currents localized at resonant surfaces. The ideal MHD models assume them so the 
opening of magnetic islands is prevented 15| , while the resistive MHD models predict them self-consistently. 
There are however other ways of plasma reaction to the perturbation, namely coupling to MHD modes which 
can provide amplification of the applied perturbation 3, 1|| 31 j, especially important at high /3 [32| . and 
those are not represented in our approach. 
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